THE SPRT FOR A BROWNIAN MOTION 


1. MAIN FORMULAS 
1.1. The general case. For a < 0 < 6 define the functions 
a=. Maes eh 
—hb — e—ha’ 

2 
T(h) = 7 (OP th) +a(1— P(h))) 

2 ae~" — be“ "4 + (b— a) 

er enhb — e—ha : 


To evaluate P, T for h = 0 we may use L’Ho6pital’s rule. We find 


Theorem 1.1. Consider a Brownian motion with drift u and infinitesimal variance 
a7. Assume we perform an SPRT test for HO: 1 = po versus H1: p = py with 
error probalities a, 8. Define 


h _ 26 — (Ho + 1) 
—————— 
M1 — Ho 
_ H1— Fo 
w= —_—_. 
oO 


Then the pass probability p and the expected duration T of the test are respectively 
given by 


p=P(h,), 
T(hy) 
T=, 


where as usual 


a=t0g (722), 
b= log (—*). 


Remark 1.2. Note that h, represents simply the relative position of w with re- 
spect to the interval [119,41] with h, = —1 corresponding to pp = po and h, = 1 
corresponding to fs = [W1. 

On the other hand w represents a universal measure for the width of the interval 
[t0, #41]. This may be motivated by the fact that the expected distance covered by 
the diffusion component of a Brownian motion is \/to. We find that, all else being 
equal, the expected duration of the test is proportional to o?. 
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1.2. The case of a symmetric test. Assume now a = £@. In that case 


‘eo 
b=—a= log ( =) 
a 


E.g. in the iconic case a = 8 = 0.05 we find 
b = —a = log(19) = 2.94444... 


In this case the universal functions P(h) and T(h) simplify as follows 


1 
PU) = Tam 
2b1—e"% 
= Ahiltenhe 


Pass probability 
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FIGURE 1. Graphs of the universal functions P and T for a = 8 = 0.05. 


2. EXPLANATION 


Consider a Brownian random walk, starting in 0 with bounds a < 0 < 8, drift 7 
and infinitesimal variance o?. Put 


h:= af. 
The probability of touching the bound 6 first is 
3 = eo 
Poi= 2 — pata 


Hence the probability of touching the bound a first is 
Da ‘= 1— po. 

By the optional stopping theorem we find that the expected time before the random 
walk touches one of the bounds is 
Paat pod 1 ae—hb — be—ha + (b— a) 

in i e—hb — e—ha ; 
Now assume that we have a Browian motion with drift 4 and infinitesimal variance 
o? and we perform an SPRT for HO: pz = wo versus H1: pw = p11. If the Brownian 
motion is at position s at time t then one computes that the log likelihood ratio is 


(ua — Ho) (2s — t(Ho + H1)) 
20? 


ie 


LLR = 
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From this one computes that LLR itself follows a Brownian motion with drift and 
infinitesimal variance given by 


_ (41 — Ho) (24 — (Ho + #1) 
om 20? , 


= 2 
ain (4 M1) ; 
(on 


We obtain 


~_ ih Hi — Ho 5 
et 2 (el 


It now suffices to substitute this into the formulas for p, and T. 


